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Anomalous dynamics of an elastic membrane in an active fluid 
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Using numerical simulations, we characterized the behavior of an elastic membrane immersed in 
an active fluid. Our findings reveal a nontrivial folding and re-expansion of the membrane that 
is controlled by the interplay of its resistance to bending and the self-propulsion strength of the 
active components in solution. We show how flexible membranes tend to collapse into multi-folded 
states, whereas stiff membranes oscillates between an extended configuration and a singly folded 
state. This study provides a simple example of how to exploit the random motion of active particles 
to perform mechanical work at the micro-scale. 


Introduction - Suspensions of bacteria and synthetic ac¬ 
tive particles offer a novel approach to manipulating 
matter at the microscale. Not only do passive micro¬ 
components in an active fluid display unusual transport 
properties [THn], but active fluids can mediate a new set 
of effective interactions between passive elements [T8H2^ , 
providing an extra handle in material design. In addi¬ 
tion, active fluids have been used to power primitive mi¬ 
cromachines [25H30]. The origin of this phenomenology 
is derived from the unique pressure (or stress) gradients 
generated by active fluids [5THiD] . 

Although, much of the effort in this field has focused 
on the interaction of active fluids with rigid, passive 
microcomponents, there has been some work on the 
behavior of flexible objects in an active fluid. Both 
flexible [111132] and semi-flexible polymers[lT| confined to 
a two-dimensional active bath have been shown to ex¬ 
hibit dynamic and scaling behavior that is much richer 
than that expected for polymers in a thermal bath. Fully 
flexible polymers [H] display a non-universal Flory scal¬ 
ing behavior as well as an anomalous chain swelling, while 
semi-flexible Dolvmers[41j display a dynamic collapse and 
re-expansion when immersed in an active bath at differ¬ 
ent values of the propelling forces. These results suggest 
that the driven nature of an active fluid can be used to 
control the shape of flexible microcomponents. This can 
be thought of as a microscopic joint or clamp where the 
mechanical action (folding) of stiff fibers can be induced 
by increasing the self-propulsion of the active components 
in solution. 

Much of the phenomenological behavior observed for 
polymers is, however, specific to two-dimensional (or 
quasi-two dimensional) systems, where active particles 
can be confined with relative ease within the bends of a 
polymer, and generate significant pressure gradients ca¬ 
pable of driving its collapse. The scaling behavior and 
statistical properties of a polymer, when embedded in a 
three-dimensional active fluid, do not differ significantly 
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from those of a polymer in a thermal bath. This would 
not be the case when replacing the polymer with a flex¬ 
ible two dimensonal surface, which is the subject of this 
work. 

Here, we study the behavior of an extended elastic 
membrane, a natural generalization of a linear polymer 
chain to an intrinsically two-dimensional structure, 
suspended in a three dimensional active bath. Physical 
examples of such materials include graphite-oxide and 
graphene sheets cross-polymerized biological 

membranes [47|, and the cytoskeleton of red blood 
cells im 135]. Polymerized membranes have been stud¬ 
ied intensively in the last few decades (see [SD] and 
references therein for a review on the subject) and 
display a phenomenological behavior that is far richer 
in complexity than that observed in linear polymers. 
Using numerical simulations, we explore the mechanical 
properties and conformational behavior of an elastic 
membrane immersed in an active fluid for different 
strengths of the bath activity and for different values 
of its bending rigidity. Our findings reveal a nontrivial 
folding and re-expansion of the membrane that is 
controlled by the interplay of its resistance to bending 
and the self-propulsion strength of the active compo¬ 
nents in solution. We show how flexible membranes 
tend to collapse into multi-folded states, whereas stiff 
membranes oscillates between an extended configuration 
and a singly folded state. 


A BCD 

FIG. 1. (Color online) Reference snapshots showing the 
various conformations taken by the membrane as the bath 
activity and bending rigidity is varied. (A) flat/extended, 
(B) bent, (C) single folded, and (D) multi-folded. 
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Methods - The elastic membrane is modeled using a stan¬ 
dard triangulated mesh with hexagonal symmetry m- 
The mesh is composed of nodes which are arranged 
in a circular geometry with diameter d. Self-avoidance 
of the surface is imposed by placing a purely repulsive 
particle of diameter a and mass m at each node of the 
mesh. Every particle on the surface is bonded to its first 
neighbor via a harmonic potential 

Ustretching — ^o) (1) 

where Kg is the spring constant, r the distance between 
two neighboring particles, and ro = the equilib¬ 

rium bond length between two particles. The bending 
rigidity of the membrane is implemented using a dihe¬ 
dral potential between adjacent triangles of the mesh; 

Ubending — COS 0) (2) 


with a range of action extending up to . Here 

Tij is the center-to-center distance between any two par¬ 
ticles i and j, and e is the interaction energy. 

Using the numerical package LAMMPS [S2], all 
simulations were carried out in a periodic box of di¬ 
mension L = 100 with T = e = m = a = T = l and 
7 = 10r“^(here t is the dimensionless time). The drag 
coefficient 7 was chosen to be sufficiently large such that 
the motion of the particles is effectively overdamped. 
The number of active particles considered in our simula¬ 
tions was N = 10, 000 giving an active particle density 
of p = N/V = 0.01. While, the number of nodes in 
the membrane was set to Nm = 1700 resulting in a 
membrane with diameter d ~ 46(7. For each simulation, 
the membrane was initialized in a flat configuration and 
the simulation was run for a minimum of 1 x 10 ®r time 
steps. All quantities in this investigation are given in 
reduced Lennard-Jones units and for convenience, we 
refer to the activity of the bath in terms of the dimen¬ 
sionless Peclet number Pe = \Fa\cT / (^D) = \Fa\iJ/{ksT). 


where (j) is the dihedral angle between opposite vertices of 
any two triangles sharing an edge and Kb is the bending 
constant. 

Each active bath particle is a sphere with mass m, di¬ 
ameter fj, and undergoes overdamped Langevin dynamics 
at a constant temperature T. Self-propulsion is intro¬ 
duced through a directional force of constant magnitude 
|Fa| and is directed along a predefined orientation vector 
n which passes through the origin of each particle and 
connects its poles. The equations of motion of an individ¬ 
ual particle are given by the coupled Langevin equations 


mir = —"fr — drV+\Fa\n+,/WDm ( 3 ) 
h = \/^r^R{t) X n (4) 

where 7 is the drag coefficient, V the interparticle poten¬ 
tial, and D and Dr are the translational and rotational 
diffusion constants respectively, satisfying the relation 
Dr = {3D)/a^. The translational diffusion constant D 
is related to the temperature T via the Stokes-Einstein 
relation D = ksT/^. The typical solvent induced 
Gaussian white noise terms for both the translational 
and rotational motion are characterized by = 0 

and = 6ijS{t - t') and {^Ri{t)) = 0 and 

= 5ij5{t — t'), respectively. Each node of 
the membrane also undergoes over-damped Langevin dy¬ 
namics at a constant temperature T where the equations 
of motion are given by Eqs. (3) and (4) while letting 
jUal = 0. The interactions between any two particle in 
the systems (membrane nodes or active components) are 
purely repulsive and are given by the Weeks-Chandler- 
Andersen (WCA) potential 


Results - To characterize and quantitatively differentiate 
between the different conformations of the membrane, 
we employ two well established shape parameters: the 
membrane asphericity A and its radius of gyration Rg. 
Following Rudnick et al. [33], we define asphericity as 
the rotational invariant 


A = 


E,>,(a.-a,)2 

(E.A.)^ 


( 6 ) 


where Ai, Aj are the fth and jth eigenvalues of the iner¬ 
tia tensor, and i,j G (1,2,3). A fully symmetric object 
such as a sphere (i.e. Ai = A 2 = A 3 ), will have a value 
of asphericity equal to zero. At the opposite extreme, 
a thin rod (i.e. Ai > 0 , A 2 = A 3 = 0 ) will have an as¬ 
phericity equal to one. We are predominately concerned 
with a planar circular membrane whose conformation os¬ 
cillates between an extended state {A « 0.25), a singly 
folded state {A « 0.5), and an ensemble of multi-folded 
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FIG. 2. (Color online) The time average of the radius of 
gyration {Rg) (left) and asphercity (A) (right) as a function 
of Pe for membranes of different bending rigidities kj,. 
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FIG. 3. (Color online) Probability distribution of the radius of gyration P{Rg) for membranes of increasingly large Peclet 
numbers. Different panels show the results for different bending rigidities of the membrane as indicated. The letters in 
parenthesis indicate the corresponding configurations as shown in Fig. 


crumpled states (that are the most symmetric) which are 
characterized by the smallest values of A. Representative 
snapshots of these configurations are shown in Fig. 
The second configurational parameter, the radius of gy¬ 
ration of the membrane, is simply the trace of the inertia 
tensor given by 

= (7) 

i 

We begin by considering the time average of the shape 
descriptors (A) and (Rg) for different values of the Peclet 
number. In the limit of low bath activity (Pe —> 0), the 
membrane is found in a characteristically extended state 
with {A) « 0.22 — 0.25 and the radius of gyration at a 
maximum {Rg) « 16 — 17cr, with the larger values as¬ 
sociated with the stiffer membranes. Notice that, unlike 
its one-dimensional counterpart (the polymer), an elastic 
membrane exhibits an overall extended state even in the 
limit for Kb = 0, which makes the dependence of Rg with 
Pe rather different than what was observed for fully flex¬ 
ible polymers confined to two dimensions [1^1112]. In fact, 


as shown in Fig. as the bath activity is increased, Rg 
systematically decreases to smaller values until a plateau 
is reached. Stiffer membranes require larger Pe before Rg 
begins to decrease and the value it decays to increases 
with the strength of the active force. The behavior of 
the average asphericity as a function of Pe, suggests a 
more complex structural landscape. While the aspheric¬ 
ity for the fully flexible membrane is relatively well be¬ 
haved where it undergoes an initial decreases followed 
by a gradual increase for larger Pe, the typical confor¬ 
mations acquired by the rigid membranes are strongly 
dependent on the value of the bending rigidity and Pe, 
as illustrated by the sharp changes in (A) particularly for 
intermediate active forces. 

To gain additional insight into the shape of the mem¬ 
brane as a function of the active forces, we analyze the 
underlying probability distributions of the asphericity 
P(A) and radius of gyration P{Rg). The results are given 
in Fig.j^and Fig. respectively. For small active forces, 
as discussed above, the membrane is nearly flat and as ex¬ 
pected the distributions for P{A) and P{Rg) are sharply 





































4 


Cl, 



0.0 0.1 0.2 0.3 0.4 0.5 

A 



0.0 0.1 0.2 0.3 0.4 0.5 

A 




A A 


^ Pe—0 

— Pe=30 

— Pe=60 

— Pe^lO 

— Pe^40 

o 

II 

1 

— Pe=20 

— Pe^50 

1 

re 

II 

oo 

o 


FIG. 4. (Color online) Probability distribution of the asphericity P(A) for membranes of increasingly large Peclet numbers. 
Different panels refer to different bending rigidities of the membrane as indicated. The letters in parenthesis indicate the 
corresponding configurations as shown in Fig. 


peaked around the corresponding values of Eg and A. 
In this low Pe limit, the distributions become increas¬ 
ingly sharper as the bending rigidity is increased. For 
larger values of Pe, both P{Rg) and P{A) broaden and 
shift towards smaller values for flexible membranes, indi¬ 
cating that the membrane is on average more compact, 
but it can also access a variety of conformations across 
the spectrum of possible shapes. More interestingly, for 
large bending rigidities distinct multiple peaks appear 
at specific values of A, which is a clear indications that 
the membrane breathes dynamically between a restricted 
number of partially stable conformations. Specifically, at 
large bending rigidities and intermediate values of Pe, 
the membrane mainly interconverts between a bent con¬ 
figuration (Fig. [^) characterized by >1 « 0.2, and a 
configuration displaying a single fold along the center of 
the membrane (Fig. A « 0.48. Visualization of the 
membrane trajectory (provided in the supporting mate¬ 
rials) reveals a continuous folding and unfolding of the 
membrane over time. For more flexible membranes, the 
dynamic behavior is similar, however, the most compact 


shape the membrane obtains is a multifolded configu¬ 
ration as shown in Fig. EP- The larger degree of flexi¬ 
bility allows for a broader range of shape deformations 
under the forces generated by the active particles. This 
is clearly observed in Fig. where we plot the joint prob¬ 
ability distribution P{A,Rg) for a fixed Peclet number 
Pe = 50 and various bending rigidities. A flexible mem¬ 
branes explores a wide range of configurations with dif¬ 
ferent Rg and A, but as the stiffness of the membrane 
is increases the explorable landscape of shape conforma¬ 
tions becomes more and more localized. 

All of the conformations active particles induce on the 
membrane surface are the result of their tendency to lo¬ 
calize near regions with large/positive curvature. The 
reasoning behind this phenomenon has been discussed in 
detail here [H1311 [3S1 m]. A rigid membrane is nearly 
flat in the absence of active particles. When active par¬ 
ticles are introduced, they tend to accumulate on the 
surface of the membrane. An initial bend of the mem¬ 
brane occurs once a sufficiently large asymmetry in the 
number of active particles develops on either side of its 
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FIG. 5. (Color online) Joint probability distribution of the 
membrane shape P{A, Rg) for a fixed Peclet number Pe = 50 
and various bending rigidities. 

area. This leads to an instability because the region of 
the membrane with higher/positive curvature will fur¬ 
ther stabilize the active particles residing on it, while the 
ones on the outer surface will be more easily escape. The 
net result is an increasing density gradient of particles 
on the surface. Similarly to the mechanism driving the 
formation of hairpins in polymers mnn], particles ac¬ 
cumulating within the creased region of the membrane 
act as a dynamic fulcrum by which the sides of the mem¬ 
brane can pivot and fold onto one another. Because no 
direct attractive interactions are present in our system, 
any folded or multi-folded configuration of the membrane 
will eventually unfold, resulting in repeating events of 
folding/unfolding transition, the rate of which increases 
with membrane flexibility. 

Two limits are worthy of a more detailed discussion. 
The first case concerns the large Peclet number limit. In 
this limit, the probability distribution P{A) associated 
to rigid membranes become comparable with that of 
their fully flexible counterparts. This is because larger 
active forces are able to fold the membrane simultane¬ 
ously at different locations, leading to the formation of 
more compact multi-folds and to an overall softening of 
the bending rigidity. Thus at large Pe the dynamical 
behavior of a rigid membrane is nearly analogous to 
that of a flexible membrane. The second limit concerns 
the case in which the membrane is sufficiently rigid that 
the active forces are not strong enough to drive the 
formation of a single fold. Here, the membrane adopts 
partially bent configurations (Fig m , and the surface 
acts as a sail able to trap a significant number of active 
particles preferentially on the side of the membrane with 
positive curvature. The net result is a active motion of 


the surface that is modulated by the propelling forces. 
The mechanism behind how an active fluid can induce 
an active transport on curved micro-components as been 
discussed in our previous work [^. 

Conclusions- In this paper we studied the behavior of 
a deformable elastic membrane in the presence of a low 
density suspension of active particles as a function of the 
mechanical parameters of the membrane and the strength 
of the bath activity. We find that as soon as the collective 
strength of the active forces becomes sufficiently large 
compared to the bending energy of the membrane a re¬ 
peating sequence of folding and re-expansion transitions 
of the surface takes place. While flexible membranes tend 
to collapse into multi-folded states, stiff membranes os¬ 
cillates between an extended configuration and a singly 
folded state. Interestingly, in the large activity-limit, the 
behavior of rigid membranes resembles that of fully flex¬ 
ible ones, indicating that strong active forces can soften 
the modes of deformation of the membrane. Further¬ 
more, we find that bent or partially folded membrane 
configurations act as sails or nets capable of trapping ac¬ 
tive particles on the positive curvatures of their surface 
and become actively transported through the medium. 

This work, is complementary to our previous study on 
rigid filaments [5], and suggests ways of exploiting the 
random motion of active particles to perform mechanical 
work at the micro-scale. The system we considered 
here consisted of a simple isotropic membrane in a 
bath of spherical particles; it is feasible to imagine that 
further control over the mechanical response of flexible 
micro-components such as the ones discussed here and 
in our previous work can be achieved by introducing 
judicious amount of anisotropy either on the elastic 
properties of the surface or on the interactions between 
active particles and surface. 
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